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1. Introduction 



Eigenvalues of random matrices have a rich mathematical structure and are a 
source of interesting distributions and processes. These distributions are natural 
statistical models in many problems in quantum physics, They occur for 

example, at least conjecturally, in the statistics of spectra of quantized models 
whose classical dynamics is chaotic, 0]. Random matrix statistics is also seen in 
the statistics of zeros of L-functions in number theory, [231 ■ 

In recent years we have seen a new development where probability distributions 
from random matrix theory appear as limit laws in models of a statistical mechanical 
nature, namely in certain random growth and random tiling problems. This came 
as a surprise and has added a new side to random matrix theory. It clearly shows 
that the limit probability measures coming out of random matrix theory are natural 
limit probability distributions. 

Only very special models, which are exactly solvable in a certain sense, can be 
analyzed. In these notes we will survey two models, random domino tilings of 
the Aztec diamond and a one-dimensional local random growth model, the corner 
growth model. We will also discuss relations between these two models. Underlying 
the exact solvability of the models is the fact that they can be mapped to families 
of non-intersecting paths and that these in turn lead to determinantal point pro- 
cesses. Point processes with determinantal correlation functions have emerged as 
an interesting class of point processes, with a rich structure and many interesting 
examples, |33| . 

2. Point processes 

2.1. General theory. We will need some general facts about point processes, 
but we will only survey those aspects that will be directly relevant for the present 
exposition, see jHJ ■ Let A be a complete separable metric space and let A/"(A) denote 
the space of all counting measures £ on A for which £(£?) < oo for every bounded 
B C A. We say that £ is boundedly finite. A counting measure is a measure £ whose 
values on bounded Borel sets in A is a non- negative integer. Typically A will be R, 
Z, some subset of these or the disjoint union of several copies of R or Z. We can 
define a er-algcbra T on A/"(A) by taking the smallest er-algcbra for which A — > £(A) 
is measurable for all Borel sets A in A. 

If B is a bounded set £(£?) is finite and we can write 

€(-B) 

(2.i) els = 
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for some xi, . . . ,x^m) G A. Note that we can have Xi — Xj for i ^ j, i.e. a 
multiple point. We say that £ is simple if £,({x}) < 1 for all x <E A. The counting 
measure £ can be thought of as giving a point or particle configuration in A. A 
point process on A is a probability measure P on J\f (A) . The point process is simple 
if P(£ simple) = 1. 

If the function <j> : A — » C has support in the bounded Borel set B we write 

(2.2) n(i+^))= na+^j)), 

i i—1 

where x\, . . . , x$(b) arc defined by 12.1(1 . If £(B) = 0, then the right hand side of 
(12.21) is = 1 by definition. Note that if \4>(x)\ < 1 for all x <E A we have 

(2.3) + fife)) = exp( / log(l + cj>{x))dt{x)). 

i J A 

A natural way to investigate a point process is to consider expectations of prod- 
ucts of the form 12.2fl . If we take for instance <j> = exp(— ip) — l,ip > with bounded 
support, we get the so called Laplace functional. We can write 

(2.4) n( 1 + ^)) = E^ E <fe)---#sO, 

where the sum is over all n-tuples of distinct points in the process and we include all 
permutations of the n points, which we compensate for by dividing by n!. We want 
to include all permutations since there is no ordering of the points in the process. If 
we have a multiple point of multiplicity k it should be counted as k distinct points 
occupying the same position. The n — term in 1(2. 4JI is = 1 by definition. Since <f> 
has bounded support, the second sum in 12.4(1 is actually finite almost surely. We 
can construct a new point process 2„ in A™ for each n > 1, by setting 

(2-5) s « = E fy*<i >■■■>*<»)> 

i.e. each n-tuple of points in the original point process £, including all permutations 
of the points, gives rise to a point in the new process. 
We define a measure M n on A™ by setting 

(2.6) M n (A) = E[E n (A)], 

for each bounded Borel set A C A™, i.e. M n (A) is the expected number of n-tuples 
of distinct points that fall in A. Here we assume that the process is such that all 
the M n , n > 1, are well-defined, M n (A) < oo for bounded A. The measure M n 
is an intensity measure for n-tuples of distinct points in the original process. The 
formula ((2.411 can now be written 

oo -. r, n 

(2.7) JJ(i + 0^)) = £ - / \{^x 3 )~ n {d n x). 



Assume that 



(2.8) > im^ Mn {B n ) < oo, 



n=0 
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where the bounded set B contains the support of 4>. Since, 



eA 1 



3=1 



< 



E 

n=0 



-M n (B n ) 



it follows from Fubini's theorem that 

(2.9) E[JJ(l + ^))] = £ri / n^)M n (d"ar). 

i n=0 ' ^ A " 3=1 

Consider the case when is a simple function 

m 

(2.10) 0(s)=5^ Oj ^(a:) 

3=1 

with Ai, . . . , A TO disjoint, measurable subsets of a bounded set B. Note that since 
the Aj's arc disjoint we have 



l + t<f>(x) = H(l+ta,) XA ^ x) 

3 = 1 



where |i| < 1, and hence 
(2.11) 



i j = l 

Set l/n!=0 if n < 0. Then, by the binomial theorem, 



3 = 1 



(2.12) 



oo m 



ETT ( ta o) nj TT 
11 n<! 11 (MA*) - m)\ 

ni,...,n m =0 j=l 3 3=1 • 7; ■ 3 ' 

°° t n 

E E 



m m 

n«?n 



. n-i . . . rim ) J (£ (A, ) — 

n=0 niH hn m =n v 7 3=1 3=1 J; 

If t, oi, . . . , a m are all positive, it follows from Fubini's theorem that 

Efjja+t^))] 



^ = E^ E 

n=0 niH \-n m —n 

On the other hand, by (|2.9|l . 



n\ . . . n Tl 



II ": - 

3 = 1 



n 



L(^)-n3-)! 



oo n „ n I m 

*ina +*(*«))] =e£t/ a n e 

i n=0 ' JA " fe=l \j=l 



(2.14) 



00 m 

E E 

n=0 ' niH h«m=n 



n 

rix ... n„ 



ajXA } (xk)j M n (d n x) 

) 'IIL 
3=1 
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Hence, for any bounded, disjoint Borel sets Ai, ... , A m in A, and m, 1 < i < to, 
such that 1 < rii < n and m + • • • + m m = rt, 



(2.15) M n (A? x • • • x = E 



n 



(e(^))! 



i - rii)\ 



This can be used as an alternative definition of the measure M n . If X is a random 
variable, E(X k ) is the fc:th moment of X, and E(X(X - 1) . . . (X - k + 1)) is 
the fc'th factorial moment of X. For this reason M n is called the factorial moment 
measure since, by Q3.8fl . can be defined using joint factorial moments. 

In many cases there is a natural reference measure A on A like Lebesgue measure 
on K or the standard counting measure on Z. We can then ask if the factorial 
moment measure M n has a density with respect to A™ on A n . 

Definition 2.1. If M n is absolutely continuous with respect to A" on A™, i.e. 
(2.16) M n (Ai, . . . , A n ) = / p n (xi,...,x n )d\(xi)...d\(x n ) 

J A±x ■ ■ ■ x A n 

for all Borel sets Ai in A, we call p n (x\, . . . ,x n ) the n'th correlation function or 
correlation density. They are also called product densities. 

We will be dealing with point processes for which all correlation functions ex- 
ist. In many cases if we are given the correlation functions (p n )n>i the process is 
uniquely determined. As can be guessed from above, the uniqueness problem is 
closely related to the classical moment problem. 

In the case of a simple point process on M we can get some intuition for the 
correlation functions as follows. Let Ai — [yi,yi + Ayi], 1 < i < n, be disjoint 
intervals. If the Ayi are small we expect there to be either one or no particle in 
each Ai. Hence, typically, the product £,(Ai) . . .£(A n ) is 1 if there is exactly one 
particle in each Ai and otherwise. From (|2.16|l we then expect 

P[one particle in each [y t ,y l + Ay t ], 1 < i < n] 

2.17) p n {y u ...,y n = hm . 

Ay^o Ay 1 . . . Ay n 

Note that p n (yi, . . . ,y±) is not a probability density. The function pi(y) is the 
density of particles at y, but since we have many particles the event of finding a 
particle at y\ and the event of finding a particle at j/2 are not disjoint even if yi ^= y 2 . 
We should think of p n (j/i, . . . , y%) as particle densities not probability densities. It 
follows from the argument above that if we have a simple point process on Z (or 
some other countable or finite set), then p n (jji, ■ ■ ■ , y n ) is exactly the probability of 
finding particles at y\, . . . , y n . 

The next proposition follows from H2.8J) and Ij2.9|l . The condition Ij2.18|l below 
implies ()2.8J) and we get Ij2.9|l . which is exactly (|3.12|l by the definition of the 
correlation functions. 

Proposition 2.2. Consider a point process all of whose correlation functions exist. 
Let <j) be a complex-valued, bounded, measurable function with bounded support. 
Assume that the support of <j> is contained in the bounded, measureable set B and 
that 



(2.18) ^^ioo J p n (xi, . . . , x n )d n \{x) < oo. 



71 = 



B" 
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M„K'x...x^)=E 



Then, 

(2.19) E [TT(l + 0(^))] = V- / Y[<f>(xj)p n (xi, . . . ,x n )d n \{x). 

TL\ /An 

Here the product in the expectation in the left hand side is defined by 12.2]) . 

We can think of the left hand side of (|2.19|l as a generating function for the cor- 
relation functions. Below we will see that l|2.19[l is useful for computing interesting 
probabilities. The condition (|2.18(l is not intended to be optimal but it will suffice 
for our purposes. 

We also have a kind of converse of proposition 12. 21 

Proposition 2.3. Let (u n ) n >i be a sequence of measurable functions u n : A" — > R. 
Assume that for any simple, measurable function (j) with bounded support, our point 
process satisfies 

oo „ n 

(2.20) E[JJ(l + 0(ari))]=£- / ]J ^K^r, . . . ,x n )dX n (x), 

i n=0 ' ^ A ™ j = l 

with a convergent right hand side. Then all the correlation functions p n , n > 1, 
exist and p n — u n . 

Proof. Arguing as above in l|2.13|) and (|2.14|) we see that 

(2.21) = / u N {x 1 ,...,x n )d n X(x). 

J A™ 1 x---xA^ 

This proves the proposition by the definition of the correlation functions and l|2.15|l . 

□ 

Proposition 12.21 is useful when we want to compute gap probabilities, i.e. the 
probability that there is no particle in a certain set. If B is a bounded, measurable 
set and l|2.18|l holds with (f> — ~xb, then 
(2.22) 

P[no particle in B] = E (Q(l - xsfo))] = Y,^T / Pn{xi, . . . ,x n )d n X(x) 

j 71=0 U ' •* ^ 

Below we will be interested in processes on R, or a subset of R, which have a last 
or rightmost particle. Consider a point process £ on R. If there is a t such that 
£(<, oo) < oo, we say that £ has a last particle. This will then be true for all t, since 
£(A) < oo for any bounded set. If x\ < • • • < x n ^ are the finitely many particles in 
(t, oo), we define x max (£) = x n r^\, the position of the last particle. The distribution 
function P[x max (£) < t] is called the last particle distribution. If E[£(i, oo)] < oo 
for some t € R, then £ has a last particle almost surely. 

Proposition 2.4. Consider a point process £ on R or a subset of R, all whose 
correlation functions exist, and assume that 

°° 1 f 

(2.23) V - / )d n X(x) < oo 
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for each iGl, Then the process £ has a last particle and 

(2.24) P[x max (0 < t] = V ^-f- / z n )d"A(z). 



n=0 



Proof. It follows from the n = 1 term in l|2.23[l that E 00)] < 00 and hence the 
process has a last particle almost surely. Take t < s. Proposition 12 . 21 implies that 



[no particle in (t,s)] = : — / 



p n {xi,. . .,X n )d n X(x). 

'(*,*)" 

We see from (I2.23|) and the dominated convergence theorem that we can let s — ► 00 
and obtain H2.24jl . □ 

Let us consider some examples of point processes. 

Example 2.5. A classical and basic example of a point process is the Poisson 
process on R with density p(x), where p is locally L 1 . Let A\, . . . , A m be disjoint, 
bounded sets in R. Then £(Ai) are independent Poisson random variables with 
parameter J p(x)dx, 1 < i < m. Hence with (ft as in l|2.10|l . 



[n( i +^))]=n E [ri( i +^) ?( ^ ) ] 




fc! 



„=o n - j=l j=l 
It follows from proposition 12 . 31 that the correlation functions are given by 

p n (x\, ...,X n ) = p(x\) . . . p(x n ), 

which reflects the independence of particles at different locations. If p(x) is inte- 

grable in [t, 00) the process has a last particle almost surely and 

(2.25) 

P[aw(0 < t\ = V y —L- / p ( Xl ) . . . p( Xn )d n \{x) = exp(- / p^cfe). 

Example 2.6. If ujv(xi, . . . ,xn) is a symmetric probability density on R^, then 
(xi, . . . , xn) — > Y]j—i S Xi maps the probability measure with density un to a finite 
point process on R. The correlation functions are given by 

Nl f 

(2.26) p n (xi, ■ ■ ■ ,x n ) = — / u N [xi,...,x N )dx n+1 ...dx N . 

(N - n)\ J R N~ n 

i.e. they are multiples of the marginal densities. This is not difficult to see using 
proposition 12.31 When point processes defined in this way are studied l|2.26|l is 
often taken as the definition of the correlation functions. 

Example 2.7. Let Hn be the space of all N x N Hermitian matrices. This space 
can be identified with M. N , since we have TV 2 independent real numbers. If p, n is a 
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probability measure on Hn and {Ai(M), . . . , Aat(M)} denotes the set of eigenvalues 
of M E Hn, then 

N 

(2.27) H N 3 M^J2S X . {M) 

j'=i 

maps /ijy to a finite point process on R. 
If dAf is Lebesgue measure ohHn, then 

(2.28) dfi N (M) = — e - TTM2 dM 

Z N 

is a Gaussian probability measure on Hn called the GUE [Gaussian Unitary En- 
semble). It can be shown, |26| . that for any symmetric, continuous function on 
with compact support 
(2.29) 

/ f{X 1 (M),...,\ N (M))dfi N {M)= / f(x 1 ,...,x N )u N (x 1 ,...,x N )d N x, 
Jh n Jr n 

where 

1 N 
(2.30) u N (x 1 ,...,x N ) = —— H i.r, r 2 \[, <' 

N ' l<i<j<N j=l 

is the induced eigenvalue measure on WL N . Hence the point process on M defined by 
(|2.3U|I has correlation functions given by l|2.26|l . 

We will show below that the correlation functions for the GUE eigenvalue process 
have a particularly nice determinantal form. This leads us to introduce so-called 
determinantal processes. 



2.2. Determinantal processes. Determinantal processes are characterized by the 
fact that their correlation functions have a certain determinantal form. 

Definition 2.8. Consider a point process £ on a complete separable metric space 
A, with reference measure X, all of whose correlation functions p n exist. If there is 
a function K : A x A — > C such that 

(2.31) p n (xi, ■ ■ .,xjf) = ^t(K{x l ,Xj)) , l :j=l 

for all xi, . . . , x n € A, n > \, then we say that £ is a determinantal point process. 
We call K the correlation kernel of the process. 

We can view the correlation kernel K as an integral kernel of an operator K on 
L 2 (A,A), 

(2.32) Kf{x)= [ K(x,y)f(y)d\(y) 

J A 

provided the right hand side is well-defined. 

Consider a determinantal process on A. Let cf> € L°°(A, A) have bounded support 
in B. Then by proposition 12. 21 
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oo „ n 

E [JJ(1 + </>(.*,))] = J2 ~ } / II ^{K{x hXj ))l j=1 d n X{x) 

j 71 = ' JA " j=1 

(2.33) =E-/ n^)^^'^))"^!^^) 



n=0 J - B ' 1 j=l 



provided 



(2.34) E^T 1 / det(K(x i ,x j ))^ 1 d n X(x)<oo. 



IV. / R „ 



The estimate (I2.34f) can usually be proved using Hadamard's inequality. The expan- 
sion 12.33fl can be taken as the definition of the Fredholm determinant 
det(J+XB^XB^)i 2 (A) = ^^{I+K(j)) L 2^ B y Here K<j> is the operator on L 2 (B) with 
kernel K(x, y)<p(y). There are other ways of defining the Fredholm determinant for 
so called trace class operators, namely det(i" + K4>))l 2 (B) = Oi(l + where {Ai} 
are all the eigenvalues of the operator K<f) on L 2 (B). If K{x, y)4>{y) defines a trace 
class operator on L 2 (B) and TrK(f> = J B K(x,x)dX(x), then theses two definitions 
agree. See |14j for more on Fredholm determinants. 

Proposition 2.9. Consider a determinantal point process £ on a subset A of K 
with a hermitian correlation kernel K(x,y), i.e. K(y,x) = K(x,y). Assume that 
K(x, y) defines a trace class operator K on L 2 (t, oo) for each t € K, and that 



(2.35) TrK = J K(x, x)dX{x) < oo. 

Then £ has a last particle almost surely and 

(2.36) P[x max (0 < t] = det(J - A r ) £ 2 (t)0o j. 

Proof. This follows from proposition 12 .41 and the above discussion provided we can 
prove l|2.34ll . Since correlation functions are non-negative and K{y 1 x) — K(x,y), 
the matrix (K{xi, Xj j) is positive definite. In that case Hadamard's inequality says 
that det(K(xi,Xj))i<i.j< n < YYj=i K{ x ji x j)- Hence 



Y.— Aet{K{x t ,x J )) l <^< n d n X{x) < ^2 — / K{x,x)dX{x) 



= exp J K(x, x)dX(x) 

by J03>- □ 

Consider for example the Airy kernel, 

poo 

(2.37) A(x,y)= / Ai (x + t)Ai (y + t)dt. 



o 



This is a Hermitian kernel and J t °° ^4(a:, :r)(£r < oo for any real t. It can be shown 
that A(x,y)x(t,oo) is the kernel of a trace class operator, and that there is a point 
process £ on K, the Airy kernel point process, which is determinantal with kernel 
A(x,y). This follows from general theory, see [331 . We have 

(2.38) F TW (t) = P[z max (£) < t] = det(7 - A) La(tj0o) . 
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The distribution function i 7 rw(i) for the last particle in the Airy kernel point 
process is a natural scaling limit of certain finite determinantal point processes. 
We call it the Tracy- Widom distribution, |36| . 

We will now look at some general ways of getting finite determinantal point 
processes. It is possible to get interesting infinite point processes by looking at 
scaling limits of these. These limiting point processes are typically of a few standard 
types, e.g. the Airy kernel point process is obtained when we scale around the last 
particle in some finite point processes on R. These infinite point processes are 
natural scaling limits and it is an interesting problem to understand how universal 
they are. In section 0] we will see the Airy kernel point process arising as a scaling 
limit of a finite point process associated with a random domino tiling of the so- 
called Aztec diamond. It also occurs as the scaling limit of GUE around the largest 
eigenvalue, see example 12. 121 

The following determinantal identity, we will call the generalized Cauchy-Binet 
identity. If we take A = {1, . . . , M}, A as counting measure on A, cj>i(k) = aik and 
ipi(k) = bik, M > N, we get the classical Cauchy-Binet identity. 

Proposition 2.10. Let (A,B,A) be a measure space, and let <j)j,ipj, 1 < i,j < N, 
be measurable functions such that (piipj is integrable for any i,j. Then, 

det ( / <fti(x)ipj(x)d\(x) 



l<i,j<N 

(2.39) = — f det(<f>i(xj))i<ij<N det(V' l (a; : ,))i< 4 j<Ard Ar A(x). 

iv ! J A n 

Proof. This is a computation, 

det I / <f>i(x)ipj(x)d\(x) ) = / det((f> i (x i )tl; j (x l )) 1 < i . : j< N d N \(x) 

\Ja / i<i,j<N JA» 



N 

1 

A™ 



4>i(xi) det(ipj(xi))i<i d < N d N X(x) 

i=l 
N 

J\ </»i(jc<r(i)) det(ip J (x a ( i) )) 1 <i J < N d N \(x) 

N 

sgn (a) Y[ i'iix^) det{ij}j(x l )) 1 < i ^< :N d N \{x) 
i=i 

I r N 

VT J AN ^2 s § n (°") XT ^(^o-(i)) det('0 J (a; 2 ))i<i J <Ard JV X(x). 



The first equality follows immediately using the definition of the determinant. In 
the third we have permuted the variables using an arbitrary permutation a € Sjy, 
and in the fourth equality we used the antisymmetry of the determinant. The 
last equality follows since the integral is independent of a. The final expression is 
exactly what we want by the definition of the determinant. □ 

Consider now the measure 

(2.40) u N {x)d N \{x) = 1 det(<fo(xj-))i<i,i<iv det(f/'i(x :) ))i<i !: ,-<Ard Ar A(x) 
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on A N , where 



(2.41) Z N =^-I det(<k(xj-))i<^-<jv det(-0 l (x i ))i< ij -< A rd Ar A(x), 

JV ! ,/ a jv 



and we assume that Zn 7^ 0. If un(x) > 0, then (|2.4U|) is a probability measure on 
A N . It follows from the generalized Cauchy-Binct identity l|2.39[l that Zm — det A, 
where A = (aj_y)i<i j<n, an d 

(2.42) ciij = J cf) i (x)ipj(x)dX(x). 

A 



Proposition 2.11. Let (A,B,X) be a measure space and let (pi^i be as in propo- 
sition \2~UK Assume that Zm given by \2.4i]) is 7^ 0. Then the matrix A defined by 
is invertible and we can define 



N 

(2.43) K N (x,y) = J2 M x )(A-%^(y). 

«,i=i 

If g € L oc {X), we have the following identity 
(2.44) 

I T\( 1+ 9( x j)) u N(x)d N X(x) = y — / TT g(x j )det(K N (x i ,Xj))i<i l j< n d n \(x), 

JA j=l n=Q JA j=l 

with un(x) given by 

Proof. That A is invertible follows from the fact that det A = Z N ^ by l|2.39(l . 
The proof of (|2.44(1 is based on the determinant expansion 

A A 

(2.45) det(/ + C) = V- Y det(C vl Ji< r . s <„, 

11—O ii,...,i n =l 

where C is an arbitrary TV x TV-matrix, and the formula 12.39fl . The identity (|2.45() 
is a consequence of multilinearity of the determinant and expansion along columns. 
It follows from (12391 and (12121 that 



A 

T[(l + g(x j ))u N (x)d N X(x) 
det(oj fc + J A (j)j{x)ij)k{x)g{x)d\{x)) 



det ^6 ]k +J2(A- 1 ) Jl jf 



det(a jk ) 

4> l (x)ip k {x)g(x)d\(x) 



(2.46) = det 5 jk / f 3 {x)h k {x)d\{x) 



A 
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where fj(x) — J2iLi(A~ 1 )ji(j)i(x) and hk(x) = g(x)t/jk(x). Using lj2.45J) we see that 
the last expression in l|2.46J) can be written 

N -, N 

,(x)hi k (x)d\(x) 



l<j,k<n 



n=0 j 1 ,...,i n =l KJA 

N 1 N If 

n=0 ' ii,...,i„ = l ' 
/ JV 

= E - / det E /i(%)^(a: fc ))d n A( a 



l<j,fe<n 

where we have used the identity Ij2.39(l in the two last equalities. Since 

N N 

^2fj(x)hj(x) =g(x) 4>j{x){A~ l )ji<i>i{y) 

we are done. □ 

Assume now that A is a complete separable metric space. If u^(x) > 0, then 
(I2.4U[1 is a probability measure on A N and the map A N 3 (x%, . . . , xn) — * Ylj=i ^ 
maps this to a point process £ on A. It follows from proposition ^ . 3l and the identity 
(|2.44|l that £ is a determinantal point process with correlation functions given by 
(|2.43|l . Although (|2.43(l gives an explicit formula for the correlation kernel it is 
rather complicated. In particular, if we want to study a scaling limit as N — > oo, 
we have to be able to find the inverse of the N x A-matrix A in a useful form. 
Sometimes it is possible to do row operations in the two determinants in H2.4()(l so 
that the matrix A becomes diagonal and hence trivial to invert. 

Example 2.12. (The orthogonal polynomial method). Consider the GUE eigen- 
value measure (|2.30|) . The density can be written as 

(2.47) u N (x u ...,x N ) = det(xr V^ 72 )^^. 

If Pj(x) is an arbitrary polynomial of degree j, j = 0, 1, ... , then by doing row 
operations in the determinant we see that 

(2.48) u N { Xl ,...,x N ) = ^ 7 ^det(p J -_ 1 (:E i )e- a!?/2 ) 2 < i;i < w . 

It now follows from proposition 12 . 1 II that the GUE eigenvalue process has determi- 
nantal correlation functions. The elements in the matrix A are given by 

(2.49) fly = / Pt-i(x)p j ^ 1 (x)e~ x2 dx. 

It is clear that it is very natural to choose pj to be the j th normalized Hcrmitc 
polynomial so that ay = 5y . The correlation kernel is then given by 

JV-l 

(2.50) K N (x,y) = P J (x)p 3 (y)e-^ 2+ y 2 ^ 2 . 
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We obtain the Airy kernel point process in the large N limit when we scale 
around the largest eigenvalue of an N x N random matrix from the GUE. More 
precisely, let A^ > X^ > ■ • • > A^ be the eigenvalues and set 

V2NXf ] - 2N 



1 ATi/3 

j > 1. Under the GUE eigenvalue measure £at = Sj=i ^ becomes a point process 
and in the limit N — > oo this process converges to the Airy kernel point process. The 
proof is based on the fact that we can investigate the scaling limit of the correlation 
kernel i|2.5U[l using asymptotics of Hermite polynomials. Also, if A max (7V) (= A^) 
denotes the largest eigenvalue then proposition 12 . 91 can be used to show that 



lim ] 



Afl/3 



(N) - 2N 



From this example, which has several generalizations, we see that orthogonal poly- 
nomial asymptotics is important in studying the asymptotics of the eigenvalues in 
some random matrix ensembles. 

2.3. Measures defined by products of several determinants. There is a use- 
ful extension of proposition ^. 1 ll to the case when the measure is given by a product 
of several determinants. Later we will see that such measures arise naturally in in- 
teresting problems. Actually both our main models will be of this type. Let X be 
a complete separable metric space with a Borel measure /x and fix m, n > 1. Fur- 
thermore, let </v,r+i :XxX->C, r — 1, . . . ,m — lbe given measurable transition 
functions, and 0o,i : Xq xI-iC, 4> m ,m+i '■ X x X m+ \ — > C given initial and final 
transition functions. Here Xo and X m= i are some given sets, which could be X or 
{1, . . . , n} for example. We will consider measures on {X n ) m of the form 

Pn,m{x)dfi(x) = 1 Wn,m(x)dfJ,(x) 

\Tl.) Zj n _ m 

(2.51) 

^ m—1 

= detiMxlx])) J] det(^ r+1 «,^ +1 ))det(0 m , m+1 (xr,< +1 ))ci At fe), 

(77..) Zj n ni 

\ / i r—1 

where x = (x\...,x m ) e (A")"\ * r = (xi, . . . ,x r n ), d^x) = UT=iU]=iM^) 
and x° £ Xq, x m+1 € X^ +1 are fixed points. Here, 

(2.52) Z n . m = ] [ w n , m (x)dii(x) 

and we assume that Z n:in ^ 0. If p n ,m(x) > we get a probability measure on (X") m . 
Set A = {1, . . . , m} x X. If we map Xj to (r, x^) € A, x gives N — mn points in A 
with exactly n points in {r} x X for each r. In this way we get a point process £ on A 
from the probability measure (|2.51|l . Let v denote counting measure on {1, ... , m}. 
We will use A = v x /x as our reference measure on A. When m = 1 we can take 
X = X m+ i = {1, . . . (/) .i(i,x) = cj>i{x) and m<m+1 (x,j) = ipj(x) to obtain 
the measure <|2.40[1 . so the present setting generalizes the one considered above. 
Our aim is to show that this more general setting also leads to a determinantal 
process. Variants of this type of setting have been developed in H3| an( i 

Ell. 
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Given two transition functions <p, if) we define their composition by <f> * ip(x, y) = 
f x (j>(x,z)ip(z,y)dfi(z). Set 

<t> r ,»{x, y) = (<f> r ,r+l * ' ' ' * <t>s-l,s){x, y), 

when r < s and <p r>s = if r > s. We assume that the transition functions are such 
that all functions 4> r ,s, 0<r<s<m + l are well-defined. This will imply that 
the integral in l|2.52|) is convergent as can be seen by expanding the determinants. 
Set A — (a.y), where 

(2.53) a id = ( j ) o, m+ i(x° l ,xf +1 ). 

Repeated use of the generalized Cauchy-Binet identity l|2.39[l gives Z n>m — det A. 
Since we assume that Z n ^ m ^ 0, we see that A is invertible. Set 

(2.54) K n . m (r,x;s,y) = K ntm (r,x; s,y) - 4> r , s {x, y), 
where r, s € {1, . . . , m\, x, y <G X and 

n 

(2.55) k n>m (r,x;s,y) = ^ <j) r , m +i (x, x" l+1 ) {A' 1 ). l0 (f> 0jS (x° , y) . 

We can now formulate the main result for measures of the form H2.51(l . 

Proposition 2.13. We use the notation above. Let g : A — > C belong to L°°(A, A) 
with support in a Borel set B C A. Let ip(r,x;s,y) = XB(x)4> r .s{x,y)g{s,y), < 
r,s < m + 1, where we omit xb(x) if r = and g(s,y) if s = m + 1. Assume 
that "0 defines a trace class operator, also denoted by ip, on L 2 (A, A) which satisfies 
Trip = f A ip(z;z)\(z). Then 



/ I| + g(r, x^))p n>m (x)dpb{x) = det(J + XBK n , m g) L 
J{x n ) m r=1 j =1 

oo _^ „ k 

(2.56) =J2~n Y[g{zj)det(K nim {z i ;z j )) 1 < i!j <kd k X(z). 

J, n * J A ^ - i 



... / a k 

It follows from proposition 12.31 that the point process £ on A is determinantal 
with correlation kernel K n ^ m (z\\ z%)^ z\, Z2 G A. 

Proof. The kernel K n m given by (|2.55|) has finite rank and by assumption ip is 
trace class. Hence XBK n ,mg is trace class and the Fredholm determinant det(7 + 
XBK niTn g) is well defined. Since also Tr (xBK n , m g) = J A (xBK n , m g)(z; z)dX(z), 
this Fredholm determinant has an expansion as given in the theorem, |14| . Write 

Z n ,m\g] = -jr-f— / JJ JJ(1 +5(^^ : ))w„, m fe)rf/z(s), 

so that Z„ im [0] = Z n ^ m = det A. Repeated use of the generalized Cauchy-Binet 
identity (|2.39() gives 

Z n ,m[g] 

= det ( / <p Q ^{x*,ti)X\{l + g{r,t r )) TT ^ r)r+1 (t r ,t r+1 )^ m>m+1 (t ro ,x7 ,+1 )d m /i(t) 



m m— 1 

/ ^O.l^i^l) II^ 1 + S{ r , t r)) II 0r,r+l(*r?*r 

J X m r=1 r=1 



/ l^Zjj^n 
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We can write 

m m 

n(l + fl(r,f r )) = l + 5^ 53 9{ri,t ri )...g(r e ,t n ) 

r=l fcl l<ri<---<r f <m 

and thus 

Zn,m[g] = det Ojj + ^2 / d l K^Osxixl M)9{nM) 

\ t=ll<n<---<rt<m^ Xt 

e-i 

(2.57) x H<t>r s ,r s+ Ats,t s +i)g(r 

s=l 

By definition 4>r,s = if r > s and hence we can remove the ordering of the r^'s in 
jjZEZjl . We obtain 



Cn m m „ 

ii)s(ri,*i) 
fc=l fclr!,...// 1 ' 



s+1, f s +l)0r £ ,m+l(^, j) 
s=l / 

= det ( Sjj + y^(A -1 )ifc / d\(r,x) / rfA(s, y)(f> ^(x° k , x)g(r, x) 



k=l 



,xj +1 ) 



(2.58) x (Er^W^)] 

where t/j*°{r, x; s, y) = 5 rtS S(x — y) and recursively 

i>* e {r,X]S,y) = / ii(r,x;u,t)ij* {e ~ 1) (u,t;s,y)d\(u,t) 



for ^ > 1. 
Set 



6(i; r, sc) = ^fc^o^a: )*;, x)y(r, af) 

fc=i 

c(r,a:;j)= / ( f>*^- x >(r, a:; s, y) J ^,m+i(y,^ m+1 )rfA(s,y) 



'A 

and let b : L 2 (A,X) — ► ^ 2 (n), c : £ 2 (n) — > L 2 (A, A) denote the corresponding 
operators. Then, by l|2.58|l . 

Z n ,m[g] 



^n,m[0] 

Now, 



det(5ij + (6c) (i, i))i<»j<n = det(7 + c6) L 2( AiA) . 



(The insertion of the xb does not change anything.) By assumption ip is a trace 
class operator and using (j) ryS = if r > s, we see that it is nilpotent, t/>* £ = if 
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l>m. Hence det(7 — ?/>) = !. Consequently, 



m— 1 



f^f = det(7 - i>) det(J + ( £ ^W)(xB*fl)) 

= det(J - V + Xs^ff) = det(J + XBKg), 
and we are done. □ 

3. Non-intersecting paths and the Aztec diamond 

3.1. Non-intersecting paths and the LGV theorem. A natural way to obtain 
measures of the form (I2.40f) and (|2.51l) is from non- intersecting paths. This is a 
consequence of the Lindstrom-Gessel-Viennot theorem in the discrete setting, 35 , 
and the Karlin-McGregor theorem in the case of non-colliding continuous Markov 
processes in one-dimension. In our applications below we will use the discrete 
setting so we will concentrate on that. 

Let Q = (V, E) be a directed acyclic graph with no multiple edges. A directed 
path 7r from a vertex u to a vertex v in Q is a sequence of vertices x±, . . . , x m in Q 
such that XiXi+i, the edges in the path, are directed edges in Q, x± — u and x m = v. 
The set of all directed paths from u to v will be denoted by LT(u, v) . If u\ , . . . , u n and 
v\,...,v n are vertices in Q, then II(u, v), u = (u\, . . . , u n ), v = (yi, . . . ,v n ), denotes 
the set of all directed paths (jri, . . . , 7r„), where 7Tj is a directed path from Uj to Vi, 
1 < i < n. We say that two directed paths intersect if they share a common vertex. 
The families of paths in II(u, v) that do not have any intersections with each other 
is denoted by II n .i. (u, v), and those that have at least one intersection by n w .;. (u, v). 
If a € S n is a permutation of {1, . . . , n} we will write v CT for (?V(i) > v a(2) >■■■■> v a(n))- 

Let w : E — > C be a given function, called the weight function, w(e) is the weight 
of the edge e in Q. The weight of a path is w(tt) = Ileeir w ( e ); i- e - the product of 
the weights over all edges in the path. The weight of paths (711, . . . , 7T„) from u to 
v is w(iri, . . . ,7r„) = w(tti) . . .w(tr n ). If S C II(u, v), then the weight of the set S 
is 

(3.1) W(S)= w(in,...,n n ). 

(7Ti,...,7r n )6S 

The total weight of all paths between two vertices u and v will be denoted by 

(3.2) 4>{u,v) = W(U(u,v)) = tt W' 

We will call <p(u,v) the transition weight from u to v. Here we are assuming that 
the sum in the right hand side of H3.2JI is convergent. We could also regard the 
weights as formal variables in some ring and (|3.3(l as an identity in that ring. 

We can now formulate the Lindstrom-Gessel-Viennot (LGV) theorem which re- 
lates weights of non-intersecting paths and determinants. 

Theorem 3.1. Let Q be a directed, acyclic graph and v = {u\, . . . ,u n ), v = 
(vi, . . . ,v n ) two n-tuples of vertices in Q such that Tl n .i. (u, v CT ) ^ only if a = id. 
Then, 



(3.3) 
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Proof. ([HED- Expand the determinant in the right hand side of Q3.3J) . By (|3.1Jl and 

det(<p(Ui,Vj)) = sgn(cr)^(ui,i; CT ( 1 ))...^(M n ,i; CT ( n )) 

n n I 

= s § n w n = e s § n n e ^ 

ctGS„ i=l oeS n i=l \ir i Gn(n i: « <T(i) ) 

= E S S n E ™(l"l , ■ ■ ■ , ^n) 

n£S„ 7ren(u,v (T ) 

sgn(cr)w(7ri,...,7r n ) + 2J E sgn.(a)w(TVi, . . . ,ir n ) 



E 

<J6S„ Trenn.i.fu.v,,) <tGS„ 7ren w . L (u,v„ 



= Sx + S 2 . 

By assumption Il n .i.(u, v CT ) = unless a — id, and hence Si = W^LTn.i^u, v)). It 
remains to show that S% = 0. 

Choose a fixed total order of the vertices, and let u> denote the first vertex in 
this order which is a point of intersection between the paths m, . . . , ir n . Let txi and 
TTj be the two paths with smallest indices which intersect w. Define a map 

(3.4) (cr,7Ti, . . . ,7T„) -> (<7',7ri, ... ,<) 

as follows. Set n' k = ir^ for k 7^ i, j, and if 

7Tj = UiXi . . . X a LUX a+ i . . . XpV a ^ 

ttj = ujyi . . . 7 Loy 1+1 . . . ysV^j) 

then 

■Ki = mxi . . . x a ujy 1+ i . . . ysVefj) 
ttj = ujyi . . . 7 ujx a+ i . . . xpv^t) . 

Also, we set a' = ao (i, j), where (i, j) denotes the transposition of i and j. Clearly, 
w(wi, . . . ,7r„) = w(tt[, ... ,<) and sgn (a) = -sgn(cr'). 

If we can show that (|3.4|l is an involution, then S2 — follows since 

S 2 = E S § n ( a ') W ( 7T '^ ■ ■ ■ ' n 'n) 

= - sgn(cr)«;(7r 1 ,...,7r T ,) = -iS , 2. 

u£S„ 7ren w . L (u,v„) 

That (|3.4|l is an involution is clear if (71^, . . . ,ir' n ) has the same first intersection 
point as (7Ti, . . . , 7r„). Since only 7^ and 7Tj were changed a new intersection point 
has to occur between them. Assume that x m is the new intersection point which 
is smallest in the total ordering. We must then have x m = X£, where one of m 
and £ lies in {1, . . . , a} and the other in {a + 1, ... , /3}, say m lies in the first set. 
But then x m . . . xi in 7^ is a cycle, which is impossible since we assumed that Q is 
acyclic. Hence (|3.4|l defines an involution. □ 

We will see below that a combination of proposition 12.111 or proposition 12.131 
with the LGV-theorem will lead us to interesting determinantal point processes in 
certain models. 
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3.2. The Aztec diamond. In this section we will discuss random domino tilings 
of a region called the Aztec diamond. The model can equivalently be thought of 
as a dimer model on a certain graph, |24| . A typical tiling of the Aztec diamond 
has the interesting feature that parts of it are completely regular, whereas the 
central part looks more or less random. In fact there is a well defined random curve 
which separates the regular regions from the irregular region, and it is this curve 
that will be our main interest. It turns out that the tiling can be described using 
non-intersecting paths in a certain directed graph, and these paths will lead to a 
description of a random tiling by a determinantal point process using the results of 
the previous section. 

The Aztec diamond, A n , of size n is the union of all lattice squares [m, m + 1] x 
[££+1], m, t G Z, that lie inside the region y\) ; \x\\ + \yi\ < n + 1}. A domino 
is a closed 1 x 2 or 2 x 1 rectangle in R 2 with corners in Z 2 , and a tiling of A n 
by dominos is a set of dominos whose interiors are disjoint and whose union is A n . 
Let T{A n ) denote the set of all domino tilings of he Aztec diamond. The basic 
coordinate system used here will be referred to as coordinate system I (CS-I). 

We can color the unit squares in the Aztec diamond in a checkerboard fashion 
so that the leftmost square in each row in the top half is white. Depending on how 
a domino covers the colored squares we can distinguish four types of dominos. A 
horizontal domino is an N-domino if its leftmost square is white, otherwise it is an 
S-domino. Similarly, a vertical domino is a W-domino if its upper square is white, 
otherwise it is and E- domino. Two dominos are adjacent if they share an edge of a 
square, and a domino is adjacent to the boundary if it shares an edge with dA n . We 
can now define four regions where the tiling has a regular brick wall pattern. The 
north polar region (NPR) is defined to be the union of those N-dominos that are 
connected with the boundary by a sequence of adjacent N-dominos, the last one of 
which is adjacent to the boundary. 

Let T £ T(A n ) be a tiling of the Aztec diamond and let v(T) denote the number 
of vertical dominos in T. We define the weight of T by letting vertical dominos 
have weight a and horizontal dominos weight 1, so that the total weight is a v ^ T \ 
If a > 0, which we assume, we get a probability measure on T(A n ) by normalizing 
this weight. When a = 1 we pick the tiling uniformly at random. 

A tiling of the Aztec diamond with dominos can be described by a family of 
non-intersecting paths. These paths can be obtained by drawing paths on the 
different types of dominos. On an N-domino we draw no path. On a W-domino 
placed so that it has corners at (0,0) and (1,2) we draw a line from (0,1/2) to 
(1,3/2), and on an E-domino in the same position we draw a line from (0,3/2) to 
(1, 1/2). Finally, on and S-domino, placed so that it has its corners at (0,0) and 
(2, 1), we draw a line from (0, 1/2) to (2, 1/2). It is straightforward to see that these 
paths form a family of non-intersecting paths from A r = (— n — 1 + r, — r + 1/2) 
to B r = (n + 1 — r, — r + 1/2), r = 1, . . . , n. The top path, from A\ to B\, can 
be viewed as a function t — > X n (t), \t\ < n, and we will call it the NPR-boundary 
process, since the north polar region is exactly the part of the domino tiling that 
lies completely above X n (i), see fig.l. 

These non-intersecting paths do not immediately, using the LGV-theorem, lead 
to a measure of the form l|2.51|l . In order to obtain a measure of this form we have 
to transform the paths. We will only outline how this is done, see [221 f° r an the 
details. Introduce a new coordinate system (CS-II) with origin at (— n, —1/2) and 
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■ I 1 ! 1 ! 1 ! 1 ! 1 ! 1 ! 




Figure 1 . An NPR-boundary process. 



axes e/j = (1, 1), fu = (— 1, 1) in CS-I, which gives the coordinate transformation 



(3.5) 



xi = x 2 - yi - n 

Ul = X2+V2- 1/2. 



In CS-II the non-intersecting paths go from Aj = (0, + to Bj = (n+l — j, —n), 
1 < j < n, and have three types of steps (1,0), (0, —1) and (1, —1), see fig. 2. We 
can view them as non- intersecting paths in an appropriate directed graph Q. The 
weight on the domino tiling can be transported to a weight on the non-intersecting 
paths by letting the steps (1,0), (0, —1) have weight a and the step (1, —1) weight 
1. Take N > n and set Aj = (0, 1 - j) and Cj = (n, -n + l-j), l<j<N, see fig. 
3. It is not so difficult to see that if m, . . . ,ttn are non-interesecting paths from 
A\, . . . , Apf to Ci, . . . , Cat, then irk has to go through Bk, 1 < k < n. Furthermore 
the paths from Bk to Ck, 1 < fe < n, and the paths from Bk to Cfe, n < k < N, 
only have steps (1,-1). Hence adding the paths from At to Ck, n < k < N, has 
no effect on the correspondence with domino tilings in the Aztec diamond or the 
weight, and we can just as well consider this extended system of paths. 

Each path iik from Ak to Ck has a first and a last vertex, which could coincide, 
on each vertical line x 2 — k, 1 < k < n. In order to get a measure of the type we 
want we have to double the vertical lines so that the first and the last vertices on 
each vertical line ends up on different vertical lines. These first and last vertices 
will form the point process we are interested in. We can also shift the paths so 
that the initial and final points, which are fixed, end up at the same height. The 
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Figure 2. CS-II and non- intersecting paths descibing the tiling. 




Figure 3. The non- intersecting paths in the graph Q. 



result can be seen in fig. 4. These non-intersecting paths which connect (0, 1 — j) 
to (2n, 1 — j), 1 < j < N, lie in a new directed graph Q' . The steps from even to 
odd columns are (1,0) with weight 1 or (1, 1) with weight a, the steps from odd to 
even columns are (1, 0) with weight 1, and we also have steps (0, —1) with weight a 
in the even columns. With this choice of weights we still have a weight preserving 
bijection with the original domino tiling of A n . The associated particles, which we 
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think of as a point process, are indicated in fig. 4. The NPR-boundary process 
corresponds to the top path in this picture. 




Figure 4. The non- intersecting paths in the graph Q' correspond- 
ing to the tiling in figure l3~^l The particles in the determinantal 
process are the circled dots. 

The paths tti, . . . , ttn just described connecting (0, to (2n, 1 < j < N , 

see fig. 4, can be thought of as being built up from 2n transition steps. We 
have points x\ , . . . , x r N on line r which connect via non-intersecting paths to points 
x 1 ^ 1 on line r + 1. Let (j) r ^ r+ i(x, y) be the transition weight to go from x 



on the line r to y on the line r + 1. It follows from the discussion above that 

a if y — x = 1 

(3-6) <t>2i,2i+i{x,y) = 



X 

1 if y — x = 
otherwise 



and 
(3.7) 



a -{y-x) if y- x < o 
otherwise 



From the LGV-theorem we see that the weight of all non-intersecting paths from 
x r 6 Z n on line r to x r+1 € Z™ on line r + 1 is 



(3.8) 



det( 



i>r,r+i{x r i,X r j ))i<i,j< 



Note that the initial and final configurations are fixed = 1 — i = xf n , 1 < i < N. 
The weight of the whole configuration of non-intersecting paths is then 



(3.9) 



det(0^ r+ i(a;[,^ +1 ))i< l!:) -<„. 



r=0 
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and normalizing we obtain a probability measure of the form (|2.51|) . We know 
that the associated point process with points (r, Xj), 1 < r < 2n — 1, 1 < j < N, 
has determinantal correlation functions with correlation kernel given by i|2.54[l . 
For reasons that will become clear below we will call this correlation kernel the 
extended Krawtchouk kernel. In order to make use of this kernel for asymptotic 
computations we have to be able to compute the inverse matrix A^ 1 in 1)2.5511 in 
some way, and obtain a more useful formula. This will be discussed in the beginning 
of the section The fluctuations of the NPR-boundary process arc described by 
the fluctuations of the particles x\, . . . , a^™ -1 , which are the last particles on each 
vertical line. If we consider a particular line, say line r, then the points x[, . . . , x r N 
form a determinantal point process and we can obtain the distribution function for 
x\, the last particle, from proposition 12.91 We will discuss the limit theorem that 
can be obtained in the section 0] 

3.3. Relations to other models. The north polar region can be investigated in 
a different way, which is related to the corner growth model that will be studied 
in section 4. This is based on the so called shuffling algorithm, [TJ]|, [T§|, which is 
an algorithm for generating a random tiling of A n where vertical tiles have weight 
a and horizontal tiles weight 1. Here is a description of the shuffling procedure 
following For a proof that it actually works, see The shuffling algorithm 
generates a random tiling of A n starting from a random tiling of A n -i. We can 
tile Ai by either two vertical dominos, with probability q = a 2 /(l + a 2 ), or two 
horizontal dominos, with probability 1 — q = 1/(1 + a 2 ). Assume now that we have 
generated a random tiling T of A n -\ according to the probability measure where 
the probability of T is proportional to a v ^ T \ Two horizontal dominos sharing a 
side of length two form a bad pair if the lower one is an N-domino and the upper 
one an S-domino, two vertical dominos sharing a side of length two are a bad pair 
if the left one is an E-domino and the right one a W-domino. Start by removing 
all bad pairs in A n _i. Next, move all remaining N-, S-, E- and W-dominos one 
step up, down, right and left respectively. After these steps what remains to fill A n 
are 2 x 2-blocks. In the final step we fill each of these 2 x 2-blocks with a vertical 
pair with probability q and a horizontal pair with probability 1 — q. This procedure 
will generate a random tiling of A n where each vertical domino has weight a and 
horizontal domino weight 1. If we draw the non- intersecting paths in a somewhat 
different way to what was done above, this shuffling algorithm can be translated 
into a certain multilayer polynuclear growth (PNG) model, see |20) . 

How does the north polar region evolve during the shuffling algorithm? It is 
clear from the description of the algorithm that it can only grow. The growth will 
be directly related to the so called corner growth model which we first define. We 
will return to this model in sectional Let A = (Ai, A2, A3, . . . ) be a partition, i.e. 
A,, i > 1 are non-negative integers, Ai > A2 > A3 > . . . and there is an £ > such 
that Ai = if i > I. The smallest such i is called the length, £(X), of the partition. 
We say that A is a partition of the integer N = Ai + A2 + . . . , written A h N. To 
the partition A we associate the following set of integer points in the first quadrant, 
the shape of A, 

(3.10) S(X) = € 7L\ ; 1 < 1 < A,, j > 1}. 
We can also define the filled-in shape of A 

(3.11) S(X) = S(X) + [-l,0] 2 
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a subset of [0, oo) 2 . The shape S(X) is one way of drawing the Young or Ferrer 
diagram associated to A. The set S(X) U (K 2 \ (0, oo) 2 ) has corners in certain 
places. These are positions where you can add a unit square to S(X) so that it 
still corresponds to a partition. Starting with the empty shape corresponding to 
A = (0, 0, 0, ... ) we grow larger shapes by adding at succesive times 1,2,... unit 
squares independently at each corner with probability 1 — q, where < q < 1 
is fixed. We call this growth model the corner growth model, JHIiCHj- Let 

<Scg(^) denote the random shape obtained at time n. 

Consider now the evolution of the NPR under the shuffling algorithm. Put 
a point in the center of each N-domino in the NPR in A n . Viewed from the 
coordinate system (CS-III) with origin at (0,n + 3/2) and axes em = (—1,-1), 
//// = (1, —1) in CST, these points form a random shape 5'a z (^) of some partition 
A. Analysis of the shuffling algorithm shows that ScG( n ) and Sa z (ji) are equal in 
law if q = a 2 /(l + a 2 ) as above. Hence results for the NPR in the Aztec diamond 
can be translated into results about the random shape in the corner growth model 
and vice-versa. 

Note that the waiting time before a corner point is added is a geometric random 
variable starting at 1. Let G*(M,N) denote the time when the point (M,N) S 
1? + is added to the square. R follows from the corner geometry that before the 
point (M,N) is added to the shape the points (M — 1,N) and (M,N — 1) must 
already have been added. Let w(i, j), i, j > 1 denote independent geometric random 
variables starting at 0, 

(3.12) F[w(i > j)=m] = (l-q)q m , 
m > 0, and define recursively 

(3.13) G(M, N) = max(G(Af - 1,N),G(M, N-l)) + w(M, N). 

It follows from the lack of memory property of the geometric distribution that 

(3.14) G*(M, N) = G(M, N) + M + N - 1. 
Consequently, we have 

(3.15) S CG (n) = {(M, N) G Z\ ; G(M, N)+M + N-l< n}. 

There is a translation of the growth rule (|3.13|) into a PNG-type model, see [3D] , 
|20| . |21j. This polynuclear growth model is not the same as the one coming from 
the shuffling procedure above. 

There is also a translation of the corner growth model into a totally asymmetric 
simple exclusion process (TASEP) in discrete time, [31], [T§|, JSj. This connection 
was used in JSj to prove the so called arctic circle theorem. The complement 
of the four polar regions in a tiling of the Aztec diamond is called the temperate 
region, and is the region where the tiling is more or less random. The arctic circle 
theorem says that, in the case of a uniform random tiling of A n , the boundary of 
the temperate region scaled down by n, converges almost surely to a circle. When 
the weight s ^ 1 we get an ellipse instead. This can be translated to a limit of 
G(M,N)/N as M,N — > oo, M/N — > const > (time constant). In section g| we 
will see how the representation of the random tiling as a determinantal process 
allows us to get precise information about the fluctuations of the arctic ellipse or, 
equivalently, of the NPR-boundary process. This then also leads to precise limit 
theorems for G(M,N). 
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The quantity G{M 1 N) also has a last-passage time interpretation, which follows 
from the recursion (|3.13|1 . An up/right path it from (1, 1) to (M, N) is a sequence of 
points (ik,jk), < k < M+N-l, such that {i ,jo) = (1, 1), (im+n-iJm+N-i) = 
(M, N) and (ik+i, jk+i)~ {ik, jk) — (1,0) or (0,1). Viewed from a coordinate system 
rotated 45° it is a simple random walk path. We have that 

(3.16) G(M, N) = max V w(i,j), 

where the maximum is over all up/right paths from (1, 1) to (M, N). 

There is an interesting limit of G(N, N) as N — ► oo if we choose the parameter 
q = a/N 2 . With this choice of q it is not hard to see, |19|. that among the w(i,j), 
1 < hJ '5= N, in each row and column there will be at most a single 1, and in 
the whole square all numbers will be < 1 with probability — * 1 as N — > oo. In 
the limit N — > oo the set {(«, j) ; 1 < i,j < N , w(i,j) > 1} scaled down by N 
will converge to a Poisson process in [0, l] 2 with intensity a. The number n of 
points in [0, l] 2 is a Po(a) random variable. Let x\ < — < x n and y\ < . . . y n 
be the x- and y-coordinates of these points. The coordinates of the points in the 
Poisson process can then be written (xj, VcrU)), j — 1, ...,fi, where a £ S n is a 
permutation of {1, . . . ,n}. The uniformity of the Poisson process implies that a 
will be a uniform random permutation from 5* n . Denote by L(a) the length of the 
longest increasing subsequence in this permutation. We say that cr(ii), . . . , o~(it) is 
an increasing subsequence in tr if i\ < ■ ■ ■ < ii and cr(«i) < • • • < cr{ii). If we look 
back at l|3.15fl we see that we should have 

(3.17) G(N, N) — ► L(a) 

in distribution as N — > oo. Hence, we may be able to use results for random tilings 
of the Aztec diamond also to study the problem of the distribution of the length of 
the longest increasing subsequence in a random permutation. 

4. ASYMPTOTICS 

4.1. Double contour integral formula for the correlation kernel. In order 
for the formulas (|2.54|l and (|2.55(l for the correlation kernel of the determinantal 
process defined by Ij2.51|l to be useful we have to find a different representation. In 
particular we need some way of computing the inverse matrix A . When A is a 
Toeplitz matrix it may be possible to do this, at least approximately, by using a 
Wiener- Hopf factorization of the symbol for A. For the models that we will consider 
A actually is a Toeplitz matrix and we will be able to find nice formulas. 

The space X in l|2.51|l will now be Z, and we will also take Xq = X m+ i. Hence 
x e (Z n ) m . Let fr(z), z = e %e , be a function in L 1 (T) with Fourier coefficients 
f r {n), neZ. Assume that the transition weights 4> r ,r+i in l|2.51|l are given by 

(4.1) (j) r>r+1 (x,y) = f r (y- x), 

< r < m, x, y G Z. Then, for r < s, 

(4-2) <p r ,s(x,y) = fr,s(V- x )> 

where 

s-l 

(4-3) frA*) = TlM*)- 

i=r 
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We see that a. tj = fo, m +i(x" l+1 - x°) = fo, m +i(i - j), since x™ +l = = 1 - j. 
The matrix A is thus a Tocplitz matrix with symbol a(z) = /o,m+i(-z); A = T n (a), 
where T n (a) = (a(i — j))i<i.j< n - A computation, |21| . shows that 

1 " 

(4.4) V k n . m (r,x;s,y)z x W -v = -/ r , n +i(-)/o,.(-) V ^(T-^a))^. 

Let T(a) = (<x(i — j))ij>i denote the infinite Toeplitz matrix with symbol a. We 
say that a S L 1 (T) has a Wiener-Hopf factorization, if it can be written a = a + a_ 
on T, where a+{z) = J2n=o a t zn > a ~( z ) = J2^=a a n z ^ n > («n)> ( a «) € ^. We can 
extend a + to \z\ < 1 and a_ to {\z\ > 1} U {oo}. We also require that a + and a_ 
have no zeros in these regions and that a has winding number zero with respect to 
the origin. Also, suppose that 



(4.5) 



for some a > 0. Then, T n (a) is invcrtible for n sufficiently large and 

(4.6) \(T-\a)) jk - (Tia+^TiaZ 1 ))^ < Cmin((n + 1 - k)~ a , (n + l-j)- a ) 

for some constant C, 1 < j, k < n. 

If f r is analytic in 1 — e < |z| < 1 + e for some e > 0, and has a Wiener-Hopf 
factorization as defined above, then a(z) = /o,m+i(z) will have a Wiener-Hopf 
factorization and (|4.5() will be satisfied. Hence we can use l|4.6[) to compute the 
inverse of the Toeplitz matrix. Combined with 14.4fl this yields 

(4.7) lim V K n<m {r,x;s,y)z x w~ v = —^—F{z,w), 

ra— >oo * — ' z — W 

for 1 — e < \w\ < 1 < \z\ < 1 + e, where 

fl\ s — 1 T — 1 77T 

(4.8) F {Z , w) =n /t - ( i) n / t + (i) n f^y n /.-^r 1 - 

t=r t=0 t=0 t=s 

Hence, 

lim K n .m(r, x; s, y) = if m (r, x; s, y) 

n — >-oo 

. t . 1 f dz f dw w y z . . 

(4-9) =TT^i\ — ~ 

\Lmy J Jr z J^ r u> z x z — w 

where 7 r is a circle with radius r and center 0, 1 — e < ri < ri < 1 + e. We obtain 

(4.10) lim K n , m (r,x; s,y) = -<f> r , s (x,y) + K m (r, x; s,y) = K m (r,x;s, y). 

n — *oo 

From (|4"2|) and it follows that 

(4.11) <j> r , s {x,y) = ±- J e^- x ^F{j\e 6 )dd. 

The integral formula (|4.9|l gives K m (r, x; s, y) if r > s. Using l|4.11|l and the residue 
theorem we see that it also gives K m (r, x; s, y) for r < s if we take 1 — e < r2 < n < 
1 + e. These integral formulas are good representations of the correlation kernel if 
we want to investigate its asymptotics. 
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4.2. Asymptotics for the Aztec diamond. In our discussion of the Aztec di- 
amond we saw that adding more paths from A k to C k for k > n had no effect on 
the probability measure since all these paths are fixed. Hence, we can take the 
limit N — > oo in the correlation kernel Kjsi jn {r,x; s,y) coming from the measure 
(|3.9|) with transition functions l|3.6[l . (|3.7() . without changing anything. We get a 
determinantal process on 2n — 1 copies of Z where the configuration is always frozen 
below the level — (n — 1). We will thus be able to use the formula l|4. 9J1 . Assume that 
< a < 1. The case a = 1 can be handled by considering the limit a — > 1 and using 
continuity. Set f2i{z) — az + 1 and /2i+i(z) = (1 — a/z)^ 1 . Then (|4.1I) holds and 
all the conditions on f r above are satisfied, /^(z) = az + 1, f%i{z) — fti+ii 2 ) ~ 1 
and f2i+i( z ) = (1 — a/z)^ 1 . Denote the generating function in (|4.8|l by F njT , iS (z, w) 
to indicate the dependence onn,r,s,0<r,s< 2n. We obtain, 

(4.12) F n> 2 r _ ei ,2 S -e 2 (^,w) - 



(l-az)"- r + ei (l + a/z) r ' 

with 6^ € {0, 1}. The Aztec diamond particle process is a determintal process with 

kernel 

(4.13) 

! C dz f dw uj y z 

K K r,n(r,X;S,y) = -<t>rA x >y) + /o -\2 / _ / 7 ^n,ry>(z, l"), 

(27Tl)^ „/ Z J W Z x Z-W 

where a < r 2 < I /a, < r\ < r 2 . Here, by (|4.11|) . 

1 f dz 

(4.14) (x,y) = 7 —— — z v*F n ^ s {z,z) 

We will call the kernel l|4.13[) the extended Krawtchouk kernel. The reason for this 
name is that it can be expressed in terms of Krawtchouk polynomials, see |22| for 
all the details. Here we will only consider the case r = s. Let pk(x;q,n) be the 
normalized Krawtchouk polynomial, i.e. it is a polynomial of degree k satisfying the 
orthogonality condition 

^2p j (x;q,n)p k (x;q,n)( n \q x (l - q) n -' x = S jk , 

< j, k < n, on {0, . . . , n}. Define the Krawtchouk kernel, 
(4.15) 



r-l 



K K ^ n ^ q (x,y) = 2_^p k (x;q,n)p k (y;q,n) 



k=0 



1/2 



It can then be shown, I22| . using the contour integral formula for pj(x;q,n) that 

(4.16) Xkr,n(2(n - r) + 1, x - r + 1; 2(» - r) + 1, y - r + 1) = K Kl ^ q (x, y), 

where q = a 2 / (1 + a 2 ). 

Note the similarity between the Krawtchouk kernel and the Hermite kernel 
(I2.5U|I . The same argument that showed that the point process defined by (|2.3U|I has 
determinantal correlation functions with kernel (|2.5(J|I shows that the probability 
measure, the Krawtchouk ensemble, 

(4.17) u n ,r(x 1 ,... > xr)=^- n (^-^) 2 n( n V(i-<zr- x ' 

l<2<j<r — 1 v J / 
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on {0, . . . , n} r defines a determinantal process on {0, . . . , n} with correlation kernel 
1)4.15(1 . Probability measures of the form ()2.30|1 and 1)4. 17j) . i.e. a Vandermonde 
determinant squared times a product of single particle weights are called orthogonal 
polynomial ensembles and go back to the early work of Gaudin and Mehta, |26|. The 
Krawtchouk ensemble is an example of a discrete orthogonal polynomial ensemble, 

m 

The NPR-boundary process is related to the top particles the 
Aztec diamond particle process. If we look first at a single line r we see that is 
the last particle in a point process on Z given by the x\ . It follows from proposition 
Eland lj4"TB|> that 

(4.18) ¥[x 2 ^ r)+1 < t - r + 1} = det(7 - Xkr,r,n,,X(t,oo))*»(Z) 

for any ( £ M. Looking at the geometry of the Aztec diamond and how the non- 
intersecting paths were defined we see that the NPR-boundary process X n (t) is 
obtained by joining the points (in CS-I) Q{j) = (2 j - x^ ax - n, x^ ax - 1/2), P(j) = 
(2j — x^i ax — n,^^ 1 — 1/2), 1 < j < n with straight lines. Hence, (|4.18|1 can be 
used to investigate X n (t). 

We see that asymptotics for the Krawtchouk kernel will give us asymptotics for 
the NPR-boundary process. The relation between the NPR-boundary process and 
the corner growth model discussed above yields 

(4.19) P[G*(M, N) < n] = F K „ l [x 2 ^~ M)+1 < n + 1 - M - N}. 
Combining (|4.18() and l|4.19|l and using l|3.14|l we find 

(4.20) P[G(M, N) < t] = det(7 - K Kl> M,t+M+N-i, q )i*at+M+i,t+M+2,... }) 
t e Z, see also [20]. 

We will not discuss the asymptotic analysis in any detail, but only give the results 
and outline the main structure of the proofs. From (|4.13() . I|4.16|l and Cauchy's 
theorem we see that the Krawtchouk kernel has the representation 

(421) (x y) 1 fdz fdw z z"~* (1 - awY(w + a)"^ 1 

1 ' K K*>r,n, q {X )y) - (2 ^ )2 J r z w z _ ww n-y (1 _ az) r {z + a) n-r+l . 

where T is given by t — > ot2 + it, t G K and 7 is a circle with radus ai centered at 
the origin, < ol\ < a-i < 1/a. 

We want to show that close to the rightmost particle the Krawtchouk kernel, 
appropriately scaled, converges to the Airy kernel given by (|2.37|l so that in the 
scaling limit close to the edge of the Krawtchouk ensemble we get the Airy kernel 
point process. The last particle in the Krawtchouk ensemble, and hence the last 
particle in the Aztec diamond point process restricted to a line, will in the limit 
fluctuate according to the Tracy- Widom distribution. In order to show this it is 
useful to have a double contour integral representation of the Airy kernel. Using 
the formula 

(4.22) Ai(x) = ^- I e iz3/3+ixz dz, 
with 77 > 0, for the Airy function in H2.37|l we obtain 

1 p p „ixz+iyw+i(z A -\-w 3 )/3 

(4.23) A{x,y) = -—2 dz dw — — . 

(2th) 2 J Imz=ri J Inlw=ri t(z + w) 
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In discussing asymptotic results for the Krawtchouk kernel we will for simplicity 
only consider the case a = 1 and the part of the NPR-boundary that lies above a 
neighbourhood of X\ = 0. 

Theorem 4.1. Set = 2- 3 / 2 (V2 + 1) and 7 = 2- 3 / 2 (V2 - 1). If r = jn, 

x = 0n + 2~ 5 / 6 n 1 / 3 £ and y = 0n + 2" 5 / 6 n 1 / 3 ?/ ; then 

(4.24) 2-^ e n 1 / 3 (V2-ir-yK Kr ^ nA/2 (x,y) - Aforj) 

as n —t 00, uniformly for £,77 in a compact set in R. Also, if X max denotes the 
rightmost particle in the Krawtchouk ensemble with q = 1/2, r = 771, i/ien 

(4.25) P[X max < 0n + 2- 5 / 6 n 1 / 3 C] -> F TW ($ 
as N — > 00 for every £. 

Proof. (Sketch) |22I - The integral (|4.21|) can be written (a = 1, which corresponds 
tog =1/2), 

(4.26) *r Kr , r ,n,i/2(z,y) = tt^I-I ^L^e^m-nf^iK 

{2-k%Y J r z w z ~ w g(w) 

where e n ^ z ^ represents the leading order behaviour of the integrand and 

/(*) = (1 - 0) logz - 7 log(l - z) - (1 - 7) log(z + 1). 

We want to apply a steepest descent argument to this integral. The saddle point 
condition f'(z) = gives (3z 2 + (27 — l)z + 1 — = 0, which has a double root if 
(7 - 1/2) 2 + 2 - = 0. This can also be written, 

(4.27) (i_ /3 _ 7 )2 + ((3 _ 7 )2 = 1/2 _ 

We see that the 0, 7 in the theorem satisfy (|4.27|l and hence we have an Airy-type 
steepest descent problem in (I4.26|) . If we carry out the argument in detail we find 
that 

2-B/«„V3( > /2_i)-* fdzfdw z cttfM _ nfM g(z) 



(27r«) 2 Jy z j 1 w z — w g(w) 

\ I' f e ixz+iyw+i(z 3 +w 3 )/3 

dz I dw- 



(27Tl) 2 7lmz=(7 Jlmw=ri l(z + W) 

as N — > 00. Together with l|4.23|l this gives (|4.24() . 

We can also use (|4.26|) to derive estimates of the Krawtchouk kernel. By propo- 
sition 12.91 the left hand side of (|4.25l) can be written as a Fredholm expansion, 
and (|4.24|l toghether with the estimates can be used to show that this Fredholm 
expansion converges to the Fredholm expansion for the TW-distribution (|2.38p . □ 



The equation 14.27fl is actually the equation for the arctic circle. If we look at 

_____ 2(71 r) I 1 

the relation 14.1811 between aWax and the position of the last particle in the 
Krawtchouk ensemble, and translate this back to our original coordinate system 
(CS-I) we obtain x\ + y\ = 1/2, which is the equation for the arctic circle. If 
we instead of a = 1 we took < a < 1, the equation for the arctic ellipse is 
x\jp + yf/q = 1, q = a 2 / (1 + a 2 ), p = 1 — q, which can be obtained similarly. 

From (|4.18() and the relation between x r ma ^ and the NPR-boundary process we 
obtain 

(4.28) P[JC n (0) < n/V2 + i^^n 1 '^] - F TW (0 
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as n — ► oo. We will consider a generalization of H4.28J) in theorem 14. HI 

From H4.26JI and H4.20|) (up to some technical details) we also obtain, in the case 
2 = 1/2, 

(4.29) lim P[G(N, N) < 2{V2 + l)N + 2 1/6 (V2 + l) 4/3 N 1/3 £] = F TW (£)- 

We will discuss a more general result in section |SJ The full result theorem 15. 31 can 
also be proved starting from l|4.20[) , (|4.21() and l|4.23() . 

Since we have the formula (|4.13() for the extended Krawtchouk kernel we should 
also be able to derive simultaneous distributions for the last particles a^ax f° r 
several lines, and hence study the convergence of the NPR-boundary process to a 
limiting stochastic process. Using (|4.16l) the limit (|4.24|) can be can be translated 
into a limit formula for KK r , n (2r + 1, x\ 2s + 1, y) when r — s. We now want to 
generalize this to the case r ^ s. The limiting correlation kernel will be the so 
called extended Airy kernel. It is defined by 

(4.30) A(r,^,rj) = e^^Ai (^ + A)Ai (, + A)dA, if r > a 

Note that A(t, £; r, rf) = A(£, rj). Using (|4.22|l it is again possible to rewrite this as 
a double contour integral. In fact, 

\ f f „ixz+iyw+i(z 3 +w 3 )/3 

(4 ' 31) ^W?L^ dZ k m dW <r-T + i(z + w) - 

where r\ > and 2r\ + r — a < in the case a > r. 

It can be shown, [3U|, [21], that there is a stochastic process, the Airy process, 
t — > A(t) with continuous sample paths almost surely such that 

(4.32) P[A(n) <&,..., A(T m ) <£ m ] =det(J-/A/) L > ({Tli ... )Tm}xR) , 

where A is the extended Airy kernel and f(jj, x) = X(£j,oc)( x ): 1 < J < m - 

We can now show the following theorem. The proof is similar to that of theorem 
14. II but somewhat more involved. It is based on l|4.13l) and (|4.31f) . see |22) . 

Theorem 4.2. Define the reseated variables ^,rj,r,a by 2r = n(l + 1/V2) + 
2 -i/6 rn 2/ 3; 2s = n(l + 1/V2) + 2- 1 /6 (7n 2/ 3; x = + 2~ 5 / 6 (£ - T 2 )n^ 3 and 

y = n/V2 + 2~ 5 / 6 (r/ - r 2 )™ 1 / 3 . Take a = 1. Then, 

(4.33) lim {V2) x ~y +2 ^e^ a - Ti ' 3+a3 ' 3 K Kr . n {2r, x; 2s, y) = A{t, £; a, rj), 

n — >oo 

uniformly for £, rj, t, a in compact sets. 

We can now give a theorem that says that the appropriately rescaled NPR- 
boundary process converges to the Airy process. 

Theorem 4.3. Let X n (t) be the NPR-boundary process and A(t) the Airy process, 
and let the weight a for vertical dominos = 1, so that we have a uniform random 
tiling of the Aztec diamond. Then, 

as ii — > oo, in the sense of convergence of finite- dimensional distributions. 
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Proof. (Sketch), [221 • The j° int distribution of the left hand side of (|4.r?4|) for 
different times r can be expressed in terms of the joint distribution of maxfc>i x r k 
with appropriate r. The fact that x r k form a determinantal process whose kernel 
is the extended Krawtchouk kernel gives that the joint distribution is a Fredholm 
determinant involving this kernel. The limit (|4.33(l and some estimates can be used 
to show that this Fredholm determinant converges to a Fredholm determinant like 
(|4.32JI involving the extended Airy kernel. This Fredholm determinant will give the 
joint distribution of the right hand side of l|4.34|) . □ 

4.3. Asymptotics for random permutations. As above we let L(a) denote the 
length of the longest increasing subsequence in a uniform random permutation from 
Sn where A is an independent Poisson(a) random variable. From 1|3.17[) we know 
that G(N, A) -> L(a) in distribution as A -> oo if q = a/A 2 . Thus by 

¥[L(a) <n]= lim P[G(A, A) < n] 

TV— <-oo 

(4.35) = lim det(7 - X Kr-iV;n+ 27V-l,a/W 2 )^ 2 ({n+7V,n+7V+l....})- 

Let ajv be given by a N /(l + a 2 N ) = a/A 2 , so that essentially on = \/ct/N. Then, 
by {OH 

A Kr ,N,n+2N-l,a/N* {x + N,y + N) 

1 f dz r dw z W y +1 {l-a N w) N (l + a N /w) N+n 
~ (27ri) 2 i 7r2 z J lri w z-w z^ 1 (l-a N z) N (l + a N /z) N + n ' 

where a at < r\ < r2 < I /on- Here we can let A — > oo and obtain 

lim K Kr Nn+2 N-x,a/N^{x + N,y + N) 

dz r dw i wy+\^ VE(w _ 1/w)+VE(z _ 1/z) 



k=0 

(4.36) 



(27Tj) 2 Z W 1 — w/z Z x+1 



= J2 J <o+k(2V^)Jy+k(^) = B a (x 7 y). 
fc=l 

The second equality follows by expanding (1 — w/z)^ 1 in a geometric series. The 
kernel B a (x,y) on £ 2 (Z) is called the discrete Bessel kernel. The limits (|4.35|l and 
d4.3(jfl and some estimates of the Krawtchouk and discrete Bessel kernels now yield 
the following theorem, |19| . 

Theorem 4.4. Let A be a Poisson(a) random variable and pick independently a 
permutation a from Sn with the uniform distribution. Denote by L(a) the length 
of the longest increasing subsequence in a . Then, 

(4.37) P[L(a) <n] = det(I - B<> ({n ,„ +li ... }) . 

The formula (|4.37(l can be used to prove a limit theorem for L(a). This theorem 
was first proved in |T| using a completely different approach. 
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Theorem 4.5. With L(a) as in the previous theorem we have 

(4.38) n L{a) al/ ^ < A = Frw® 
as a — » oo . 

Proof. (Sketch). The Bessel functions have the following asymptotics 

(4.39) a 1 / 6 J 2v ^ +?Ql/6 (2>/5) -> Ai (0 

uniformly for £ in a compact interval as a — > oo. This, together with appropriate 
estimates, gives 

/>oo 

(4.40) a 1/6 B a {2 v ^ + ^a 1/6 ,2y^+r ] a 1/6 ) -» / Ai (f + t)Ai (r? + t)d* = 77) 

j 

as a -> 00. Thus, by l(Q7Jl . 

(4.41) lim P[ L(a) ~ 2%/ " < t] = det(J - A) i2(t;0o) - F TW (t). 

a^oo Q! ' 

□ 

If £jv(c) denotes the length of the longest increasing subsequence in a uniform 
random permutation from Sn, a de-Poissonization argument, |17| . pQ, makes it 
possible to deduce 

from the previous theorem. 

5. The corner growth model 

5.1. Mapping to non-intersecting paths. In this section we will give another 
approach to the corner growth model by mapping it to non-intersecting paths in a 
different way than that related to the Aztec diamond. 

Consider a right/down path A from (0, L) to (K,0), i.e. a sequence of points in 
Z 2 , P 3 (X) =(x j ,y j ),j = 0,...,K + L, such that P (A) = (0, L), P k +l(X) = (K, 0) 
and Pj+i(A) — Pj(X) = (1,0) or (0,-1). We use the same notation, A, as for a 
partition since there is a unique associated partition, A& = max{xj ; yj = k — 1}. 
Note that Ai = K and = L. 

Let w(i,j), <E Z 2 ^, be independent geometric random variables with param- 
eter a,ibj, < ai, bj < 1, 

(5.1) P[w{i,j) = m] = (1 - a^OM,)™ 

to > 0. Also, we set w(0, j) — w(j, 0) = 0, j > 0. Define as previously, l|3.13fl . 

(5.2) G(i,j) = max(G(z - 1, j),G(i, j - l))+w(i,j), 

€ Z2_, where G(i,0) = G(0,i) = 0, i > 0. 
Given a partition, or right/down path A, set W(A) = (w(«, i))(i,j)eS(A) > where 
5(A) is the shape of A as defined previously. If A = (M, . . . , M, 0, . . . ) with £(X) = 
N, then W(X) is the M x TV-matrix {w{i, j))\<i<Ms<j<N ■ We want to map W(A) 
to a family of non-intersecting paths in a weighted graph G(X), in such a way that 
the top path gives the values (G(Pj(X)))^ L of G(i,j) along the right/down path. 
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Corresponding to l|5.1|l we define the weight of W(A) by 

(5.3) n i^r (i ' j) - 

(M)eS(A) 

We also want the mapping to the non-intersecting paths to be weight-preserving so 
that we can use the paths to study G(M, N). 

The directed graph G{X) is defined as follows. The vertices are {— L, — L + 
1, . . . , K}xZ, and the undirected edges connect (i, j), j) for i = —L,... f K—l, 
j 6 Z (horizontal edges) and (i, j), + 1) for i = —L + 1, . . . , K, j € Z (vertical 
edges). The step Pj(X)Pj + i(X) in the right/down path A is a right-step if Pj{\) = 
(i — 1, x) and P J+1 (A) = (i. x). In that case the vertical edges with first coordinate 
— L+j + 1 are directed from (— L+j + 1, k) to (— L+j + 1, fc + 1), k e Z, i.e upwards, 
and are given the weight a, (up-edges). The step Pj(A)Pj + i(A) is a down-step if 
Pj(A) = (a;,z) and Pj+i(A) = (a;,i — 1). In that case the vertical edges with first 
coordinate —L+j+1 are directed from (— L+j + 1, k) to (— L+j + 1, k— 1), fc G Z, 
i.e downwards, and are given the weight bi {down-edges). All horizontal edges are 
directed to the right. 

To a path it in G(X) from (— L, 1 — j) to (If, 1 — j), for some j > 1, we can 
associate points Qi(n) — (i — L,Xi), < i < K + L. We let Qi(ir) be the last vertex 
in the directed path with first coordinate — L + i. If the edges on the vertical line 
x = — L + i are up-edges then Xi > Xi_i, if they are down-edges, then Xi < Xi-\. 
We have xq = xl+l = 1 — j- 

We can now formulate a theorem which gives the mapping from W(A) to non- 
intersecting paths in G(X). To our knowledge this theorem has not appeared in 
its present form in the literature so we will give a proof in section [SJ For an 
investigation of the behaviour along right/down paths in the Poissonian case see 

M 

Theorem 5.1. Let A be a partition giving a right/down path {Pj(\))f+ L , L = t(\), 
K = Ai, from (0, L) to (K,0). There is a one-to-one weight preserving mapping 
from W(X) with weight to non-intersecting paths (tti,it2> ■ ■ ■) w the weighted 

directed graph G{X), where Ttj goes from (—L,l — j) to (K,l — j), j > 1. The 
path TTj consists only of horizontal edges if j > min(if, L). If Qi(iri) = (i — L, xi), 
< i < K + L, are the points associated to the top path iri, then Xi = G(Pj(A)). 

To the paths tti, 7T2, . . . we can associate a point configuration (r, Xj), — L < r < 
K, j > 1, by letting Qifa) = (i = L, x i f L ). Note that xj L = xf = 1 - j, j > 1 are 
fixed. Also, Xj = 1 — j, —L <r<K,ifj> min(K,L). Hence theorem 15 . II maps 
W(X) with probability measure l|5.1|l to a point process in {— L + 1, . . . , K — 1} x Z. 
By the general formalism presented above this will be a determinantal point process. 

Theorem 5.2. Let X be a partition, a right/down path from (L,0) to (0,K), 
L = £(X), K = X\. The probability measure 15. 1\) on W{X) can be mapped to a 
determinantal point process on {—L + 1, . . . , K — 1} x Z. Set 

(5-4) /+( Z ) = I^L , f-(z) = l 

1 — OiZ 

if the edges on x — r + 1 are up-edges with weight ai, and 
. . . . . . . . 1-6,- 
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if the edges on x = r + 1 are down- edges with weight bi, —L < r < K. The 
correlation kernel is then given by 

\ f dz f dvj nfl z 

(5.6) Kx(r,x;s,y) = -(f> riS (x,y) + 7 - / — / F x (z,w), 

{2ni) 2 ] lr ^ z J w z x z— w 

where max(&i) < r\ < T2 < min(l/aj) ; 

A'-l s-1 r-l A-l 

(5.7) f X {z,w) = n /,-(-) n /* + (-) n n /rc-r 1 . 

r s = if r > s and 

(5.8) <Mz, V) = ^f e^-^F x (e ie ,e l9 )dd 
if r < s. 

Proof. (Sketch) The transition weight to go from (r, x) to (r + l,y) in C/(A) is 

(5-9) <f> r , r+1 ( X ,y) = Q iiy ~ X 

10 it y < x 

if the vertical edges on the line r + 1 are up-edges with weight Oi, and 

if > 2/ 



(5-10) (t> r ,r+l(x,y) 



if a; < y 



if they are down-edges with weight bi. It follows from theorem 15. II and the LGV- 
theorem that the probability measure on the point configuration is given by 12.51fl . 
Hence it has determinantal correlation functions by proposition 12.131 We have 
<fi r . r +i(x, y) — f r (y — x) , where f r = /+ f~ . The matrix A in l|2.55|l is a Toeplitz 
matrix and we can use the Wiener-Hopf factorization technique in section 0] to see 
that the correlation kernel is given by (|5.6|l . compare (|4.9|l . □ 

5.2. The Schur and Plancherel measures. Take A = (N, . . . , N, 0, . . . ), £(\) 
X . Then we have only up-edges to the left of the origin and down-edges to the 
right of the origin. Furthermore G(N,N) = maxj>ii". Restrict the attention to 
the point process x®, j > 1, above the origin. This is then a determinantal point 
process with correlation kernel 

1 f dz f dww v z JX- (1 - aj/z)(l - bjw) 

(5 - n) Kn{x ' v) ~w? L 2 t y 7ri --— ii (i - aj/W )(i - h z) ■ 

Consider now the case a, = y/q, 1 < i < N, bi = y/q, 1 < i < M, bi = 0, 
M < i < N, where M < N. The kernel can then be expressed in terms of Meixner 
polynomials, |18| . and is called the Meixner kernel. If we scale appropriately around 
the last particle this kernel has the Airy kernel as its scaling limit and proceeding 
as we did for the Krawtchouk ensemble we can prove the following theorem, 18, 

Theorem 5.3. Let < q < 1 and consider G(M, N) defined by where w(i,j) 

are i.i.d. geometric random variables with parameter q. Then, for 7 > 1, 

mt G([>yN],N)-NL>h,q) . , _ , . 

hm P . r < s = Ftw{s) 

a^oo N 1 / 3 a(^,q) 
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for any s € K, where 



, v (1 + v^) 2 
LO{"f,q) =- 



1-3 



-(7, g) = g ' 1 1 /6 _ 7 g 1/6 (V7 + V^) 2/3 (i + V^) 2/3 - 

The probability measure on the points x\ > ■ ■ ■ > £jy (M = TV) is the so-called 
Schur measure, [2Hj ■ The points x\ > ■ ■ ■ > x^j can be related to a partition p, 
by x® = Hi — i + 1 so we can also think about the Schur measure as a measure on 
partitions. This is the probability measure on partitions p that we obtain if we map 
(w(i,j)), an TV x N matrix, w(i,j) as in (53) , to a pair of semi-standard Young 
tableaux with shape \i using the RSK-correspondence, [22] • We will not prove this 
here. In the present case all the vertical edges to the left of the origin and at the 
origin in the directed graph are up-edges, whereas those to the right of the origin 
are down-edges. The non-intersecting paths pass through the points pi — i + 1, 
1 < i < N, above the origin. By the LGV-theorem the Schur measure can thus be 
written 

(5.12) p s (p) = — det(0(l -j tfH -i+ 1)) det(V>(l -j,Pi~i + 1)), 

AN 

where 

< 7j( M ; u )= X! ail ■■■ a iv-u — h v - u (ai,...,a N ), 

l<ii<---<i v - u <N 

ip(u,v)= ^2 b il ...b iv _ u =h v _ u (b 1 ,...,b N ) 

l<H<—<i v ^ u <N 

Here hk(ai, . . . , ajv) is the fc'th complete symmetric polynomial. The symmetric 
polynomial 

(5.13) s^ai, . . - ,a N ) = det(h fii - i+j (ai, . . . ,a N )) 

is the Schur polynomial labelled by ji. The Schur measure can thus be written 

(5.14) Ps(m) - ^--V( a i' • • • ) «Jv) s m( & 1' ■ • • ) b N), 
which explains the name. The normalization is, |32|. 

N 1 

(5.15) Z N = y]sn(ai, ... ,ajy)an(6i } .. . , fry) = " 



It follows from above that, under the Schur measure, X^t>i ^Ai-i is a determinantal 
point process on Z with correlation kernel 1|5.11JI . 

If we restrict our attention to the case when W — (w(i,i))i<i,j<jv is a permuta- 
tion matrix, then G(N,N) — £n(o~) is exactly the length of the longest increasing 
subsequence in the permutation corresponding to W. Restricting to a permutation 
matrix means that we want to have exactly one up-step of size 1 on each line with 
up-edges and exactly one down-step of size 1 on each line with down-edges. If the 
height configuration at the origin is p,i — z + 1, i > 1, £n(g) = fix, the corresponding 
measure on partitions p is 

1 

JV! 1 



(5.16) Ppian,Jv[M] = T77[ a l ■ ■ ■ • ■ • M S M( a )v(k)> 
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where [ai . . . a^bi . . . 6jv] means that we take the coefficient of the monomial 
ai . . . a/v&i • • • &w in S/i(a)s M (6). This measure is called the Plancherel measure, 
and is also given by f^/N\, where / M is the number of standard Young tableaux 
with shape fj,. We saw above that when studying the problem of the length of the 
longest increasing subsequence in a uniform random permutation from Sn it was 
natural to let TV be a Poisson(a) random variable. We thus consider the Poissonized 
Plancherel measure, 



a N 



(5-17) P?pM= E M Ppian ' JvM 



where Ppi an ,iv[/*] = if N is not a partition of AT. 

Theorem 5.4. d^EO^ Under the map fi — > 2j>i *> ^ e Poissonized Plancherel 
measure is mapped to a determinantal point process with correlation kernel B a (x, y), 
the discrete Bessel kernel, given by j^.^6| ). 

Proof. (Sketch). Suppose that <? : Z — * C has support in [— L, oo) U Z for some 
L > 0. By fCTl - (fS~TTI) . 

OO 

i" i=i 

oo _/v iV ^ oo 

( 5 - 18 ) =Efp[«i-^'-M II T^TTrEIKi+^-i))^^)- 

7V=0 ^ i,j=l 1 3 M j=l 

Now, since the Schur measure has determinantal correlation functions we know that 



(5.19) 

oo oo _^ k 

^2Yl( 1 +9(fJ'j-j))ps(fJ')=^2-j^ E n9(zj)<tet(.KN(xi + l,Xj + l))i<i,j<k, 

M 3=1 fc=0 ' xe1 k 3=1 

where if^ is given by (|5.11|l . Inserting l|5.19|l into l|5.18|l a rather long computation, 
which we omit, gives 

oo 

E lit 1 + 9(Pi - i))PppN = det(J + S Q 5 )p (z) , 

I" 3=1 

which proves the theorem by proposition ^. 31 □ 

5.3. A discrete polynuclear growth model. Consider now the right/down path 
which is given by A = (N, N — 1, . . . , 1, 0, . . . ) and choose a, = 6, = ■ s fq. set 
G(i + 1/2, j + 1/2) = G(i,j), i,j > 0, and 

, t + x + 1 i - a; + 1 , 

(5.20) %, t) = G( , ) 

for x e Z, i > 0, |x| < t and fi,(x,t) = if \x\ > t. We see that h(x,N), x = 
—N,...,N are exactly the values of G(i,j) along the right down path A. Set 
u(x,t) — if t — x is even or if |a;| > t, and 

, , , > ,t + x + 1 t- x + 1 

(5.21) u}(x,t) = w( , ) 
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otherwise. From (15.2(1 it follows that 

(5.22) h(x, t + l) = max(h(x - 1, t),h(x, t), h(x + 1, t)) + u(x, t+1). 

We think of h(x, t) as the height above x at time t. This growth model is a dis- 
crete polynuclear growth model, |25| . It is different from the one related to the 
shuffling procedure in the Aztec diamond. It follows from theorem 15.21 that the 
height fluctuations above [ct], < c < I at time t are of order t 1 / 3 for large t and 
are described by the Tracy- Widom distribution. But we actually know the whole 
extended kernel l|5.6l) for A = (N, N — 1 , . . . ) and hence we can study x — ► h(x,t), 
\x\ < t, as a process. This was done in j^O] for the Poissonian limit of the model 
corresponding to random permutations, where we have an extended discrete Bessel 
kernel. 

If we consider the process at time 2N — 1, and look at even x, we are studying 
the process u — > G(N + u,N — u), \u\ < N. The appropriate scaling limit of the 
kernel (|5.6|l will again be the extended Airy kernel, similarly to what we got when 
we studied the NPR-boundary process. We have the following theorem, [23, which 
we will not prove. Basically we have again to investigate the convergence of the 
extended kernel l|5.6|l to the extended Airy kernel using the saddle-point method. 
We define the rescaled process ifjv(i) by 



(5.23) G(N + u,N-u)= N + dN x '*H N (- — ^dN~^ 3 u ) 



where 
(5.24) 

i - q 

We extend it to a continuous process defined for all times by linear interpolation. 

Theorem 5.5. We have H^(t) — > A(t) — t 2 , where A(t) is the Airy process, in the 
sense of convergence of finite- dimensional distributions. 

From (|3.I6I) we know that G(N, N) is a certain point-to-point last passage time. 
It is also natural to consider a point-to-line last passage time, 

(5.25) G-pi(N) = max G(N + u,N-u), 

\u\<N 

i.e. we take the maximum of over au up/right paths from (1,1) to the 

line x + y = 2N. A maximal path will be like a directed polymer with one free 
end. From (|5.23|l we see that it is natural to consider the maximum of the process 
-ff]v(t). To do so we need a stronger form of convergence then in theorem 15 .51 The 
next theorem is proved in |21| . 

Theorem 5.6. There is a continuous version of the Airy process A(t) and H^ft) — > 
A(t) — t 2 as N — > oo in the weak star topology of probability measures on C{— T, T) 
for any fixed T . 

It is proved in jjJJ that 

(5.26) P[ PU ' ^i /3 ^ q> < s] - Fi (s) 

as TV — * oo, where Fi is the largest eigenvalue, Tracy- Widom law, for the Gaussian 
Orthogonal Ensemble (GOE). This law is different from Ftw, which is often de- 
noted by F 2 , |37|- We will not give its explicit form here. If we combine this result 
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with theorem 15.61 we see that 

(5.27) F 1 (s)=F[sup(A(t)-t 2 ) < s\. 

ten 

It would be interesting to have a more direct approach to this result. 

The maximal path in the point-to-line last passage problem is not necessarily 
unique so there could be several possible endpoints on the line. Set 

(5.28) K N = inf{s; supH N (t) = supH N (t)}, 

t<s t£R 

and similarly for the limiting process 

(5.29) K = in£{s ; aup(A(t) - t 2 ) = sap{A{t) - t 2 )}. 

t<s teR 

If we could show that the process t — > A(t) — t 2 has a unique point of maximum 
almost surely it would follow that — > K, and the law of K would be the law 
of transversal fluctuations of the endpoint of a maximal path in a point-to-line 
problem. However, the above argument gives us no clue what this law could be. 
For all we know it could be Gaussian. 

5.4. Proof of theorem l5.ll The procedure is very close to the Robinson-Schensted- 
Knuth correspondence. See [2Z| for a related analysis. 

We have to define the paths 7ri,7r 2 , .... Set w^°'(i,j) — w(i,j). Assume that 
we have defined w^(ij) for some k > 0, i, j > 0. Set G^(0, j) = G (fc) (j,0) = 0, 
j > 0, and 

(5.30) G« = max(G« (i-lj), G (fe) (i, j - 1)) + 
i,j > 1, so that G (0) (i, j) = G(i,j). Also, define 

(5.31) w^frj) - min(GW( l - l,j),G^(i,j - 1)) - G«(i — l,j — 1). 

This defines w^(i,j) and G^(i,j) recursively for all k > 0, i,j > 1. 

The path n/. goes between the points (— L + j,G^ k ~ 1 '(Pj(X)) — k + 1), j = 
0, . . . , K + L and respects the direction of the edges in £7 (A). 

Claim 5.7. We have that w^(i,j) = 0if0<i<kor0<j<k. 

To prove the claim we use induction on k. The claim is true by definition if 
k = 0. Now, if wW(i,j) = for < i < k or < j < k, then G^(iJ) = for 
< i < k or < j < k. If k < k + 1 or j < k + 1, then i - 1 < k or j - 1 < k. 
Hence, - l,j) = or G^(i,j - 1) = 0, so u/ fc+1 )(i, j) = by lESTI . This 

proves claim IBT7I 

Claim 5.8. The paths 7Ti,7T2, ■ • • in £?(A) c?o noi intersect. 
We first prove that 

(5.32) min(G (fc) (i - 1, j), G (fc) (i, j - 1)) > G (fe+1) («, j) 

for ali i, j > 1, fc > 0. If = (1,1), then fiHty is clearly true since G (fe+1) (1, 1) = 

for all fc > by claim E2| Assume (A) that lf5~52|) holds for 1 < i, j < m. If 

1 = m + 1 and j = 1, then (fK~TT5f) holds because G (fc+1) (m + 1, 1) = by claim l5~7l 
Assume (B) that Ij5.32|l holds for i = m + 2, 1 < j < n, where n < m. We want to 
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prove l|5.32jl for i = m + 1, j = n. By Ij5.30jl and [|5.31jl the inequality (|5.32() with 
i = m + 1, j = n is equivalent to 

min(G (fc) (™ + 1), G (fc) (m, n)) > max(G (fe+1) (m, n), G (fc+1) (m + 1, n - 1)) 

+ min(G (fe) ( m + 1 , n — 1 ),G (k) (m,n))-G {k) (m,n-l) 

i.e. 

(5.33) G (fe) (m, n - 1) > max(G (fc+1) (m, n), G (fc+1) (m + 1, n - 1)). 

By assumption (A), G^ fc ' (m, n— 1) > G' fe+1 ) (to, n) and by assumption (B) (m, n— 
1) > G( fc+1 )(TO + l,n- 1), so fESB holds. 

By induction l|5.32() also holds for i = m + 1, 1 < j < to. A similar argument 
shows that (I5.32|) holds for 1 < i < to, j = to + 1. It remains to consider the case 
= (to + 1, to + 1). The H5.32f> is equivalent to 

(5.34) G [k) (to, to) > max(G (fe+1) ( TO, TO + 1 ),G (fc+1) (TO+ 1,to)) 

by (|5.3()jl and (|5.31jl . Now i = to, j = m + 1 in Ij5.32jl gives G (fe ^(m, to) > 

G*( fe+1 )(TO,TO+ 1), and i = m+ 1, j = m in (|5~^|l gives G^(to,to) > G( fc+1 )(TO + 
1, to). Hence l|5.34l) holds and we have shown that (|5.32(l holds for 1 < i, j < m + 1 
and hence holds for all i,j > 1 by induction. 

We can use (I5.32f) to prove claim IS~%1 Consider an up-step Pj{\) to P J+ i(A), 
Pj(A) = (i-l,x),P J+1 (X) = (i,x). The path 7T fe then goes from (-L+j.G^-^CPj (A))- 
fc + 1) to (-L + j + l 7 G^ k -^{P j+1 {X)) -k + 1) via the points (-L + j + 1,to), 
where G^^P^A)) - k + 1 < to < G^" 1 ) (P j+1 (A)) - fc + 1. Non-intersection 
between ir k and 7r fc+1 holds if G^^P^A)) - k + 1 > G^(P j+1 (X)) - k or 
G( fe-1 )(i — l,x) > G( fe )(«,x), which follows from (|5.32|) . A down-step is treated 
analogously. This proves claim l5~%l 

Claim 5.9. The path nk is horizontal if k > min(K,L). 

Assume that K < L and k > K. the path tt^ goes through the points (— L + 
j,G( fc - 1 )(P 7 (A)) -k+1). here P,(A = (x,y) where x < K. Now, G( fc " 1 )(x,j/) = 
if x or y is < k by claim l5~?l so irk goes through the points (— L + j, —k + 1) and is 
horizontal. 

Claim 5.10. Set W^{X) = (w^ (i, j))(ij)es(A) and define the weight ofW^(X) 
to be 

n M,r (fc,( - ) . 

Then, the weight of W^^ 1 ^ (A) is equal to the weight of W^ k \X) times the weight 
of TTk in G(\) for all k > 1. 

Let 1 < i < K . In the weig ht of W( k -V(\), ai enters as 

if the up-step with weight is the step from (i — 1, to) to (i, to) in the right/down 
path A. The up-step has size G^" 1 ) (i,m) — G^^^i — l,m) so we want to show 
that 

m m 

(5.35) G^" 1 ) (i, to) - G (fe - 1} (t - 1, to) + w (fe) (i, j) = ^ w (fc - 1} (t, j). 
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Now, by ffT5n)l and (fHTTT^i . 

w (k -^ (i, j) = (i, j) - max^- 1 ) (i - l,j),G^ (ij-l)) 

-w<V{i,j) = - 1,3 - 1) - min^" 1 ^ - 1, j), G^ (i, j - 1)). 

Adding these two inequalities and summing over j we obtain 

m m 

(h 3) ™ {k) (i, i)) = El G(fe_1) ( f > J) - G(fe_1) ( z - *> 3)] 

3=1 3=1 
m— 1 

- E & J) - G( ^ 1} (* - !» J')] = G(fc ~ 1} (*'' m ) - G(fc_1) (*' - 1> m - 1), 

so we have established ()5.35f) . The argument for 6^ is analogous and we have proved 
claim QUI 

To prove the theorem it remains to show that the mapping from W^A) to 
7T2, . . . ) is one-to-one and weight preserving. That the map is weight preserving 
follows by claim 1531 and repeated use of claim IOTjI We have to show that the map 
is invertible, i.e. given (7Ti, 7T2, . . . ) we can reconstruct W{\). 

If k > win(K, L), then W {k) (A) = by claimO Assume that we know W (k) (A) 
for some k > 1. We want to show that we can reconstruct W^ 1 ^ (A). Repeating 
this we eventually get W^\\) = W(X). 

Let A be a partition of TV and let A = A° > A 1 > ■ • • > \ N ~ 1 > X N = be 
a sequence of partitions such that we get X k from A fe_1 be removing one point on 
the boundary. (We have a path from A to in the Young lattice, [62].) Assume 
that we know all the values of G^ -1 ) (m,n) for (m,n) along the boundary of A J . 
From 7Tfe we know all the values of G^ k ~ x > (m,n) along the boundary of A = A. 
Let (to, n) be the point on the boundary of A J that we remove when we go to 
A J . Since A J+1 is also a partition, (to — l,n) and (m,n — 1) must also be points 
along the boundary of A J . Hence, we know G(* -1 )(m — l,n) and G^ -1 - 1 (to, n — 1). 
We can now get (to, n) from l|5.30|) . Since we know uA fc ) (to, n) we can also 

get G^ _1 - ) (to — l,n— 1) from (|5.31(l . hence we know the values of G^ -1 * 1 (to, n) 
along the boundary of X^ +1 . Proceeding in this way we succesively get the numbers 
ui( fe_1 ) (to, n) for (to, n) G S'(A). This completes the proof of theorem l5.ll 
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